Absence of Wigner Crystallization in Graphene. 
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Graphene, a single sheet of graphite, has attracted tremendous attention due to recent experiments 
which demonstrate that carriers in it are described by massless fermions with linear dispersion. In 
this note, we consider the possibility of Wigner crystallization in graphene in the absence of external 
magnetic field. We show that the ratio of potential and kinetic energy is independent of the carrier 
density, the tuning parameter that usually drives Wigner crystallization and find that for given 
material parameters (dielectric constant and Fermi velocity), Wigner crystallization is not possible. 
We comment on the how these results change in the presence of a strong external magnetic field. 

PACS numbers: 



Graphene is a single sheet of carbon atoms which can 
be extracted from graphite by micro-mechanical cleav- 
age Various theoretical and experimental studies of 
this 2-dimensional system have shown that its properties 
H H H H E Q are markedly different from those of the 
conventional 2-dimensional electron gas (2DEG) formed 
in semiconductor heterostructures. The carriers in such a 
2DEG are well described by Fermi liquid theory, and have 
quadratic dispersion at small momenta, = Ti'lP' jlm* , 
where m* is the band mass 8]. Carriers in graphene, on 
the other hand, have a band structure in which electron 
and hole bands touch at two points in the BrouUion zone. 
This, combined with the hexagonal lattice structure of 
graphene, lead to linear dispersion at small momenta, 
-Efc = hvck. Therefore the carriers behave as "massless" 
Dirac particles with characteristic velocity 

VG ~ 10^ m/s. This difference in the character of quasi- 
particles - massive quadratically dispersing vs. massless 
linearly dispersing - gives rise to a host of remarkable 
phenomena including unusual quantum Hall effects and 
TT Berry phase 

Another interesting aspect of graphene is that it is 
a better realization of two dimensional system. The 
conventional 2DEG is formed in semiconductor quan- 
tum wells where the width of the quantum well is usu- 
ally around 100 - 300 A. Various properties of such 
a 2DEG have been extensively studied using primar- 
ily transport measurements. However, it has not been 
amenable to local probes such as a scanning tunneling 
microscope (STM), because the 2DEG is buried about 
1000 A from the sample surface. In graphene, on the 
other hand, the width of the effective quantum well is 
approximately 5 — 10 A (distance between two graphene 
sheets in graphite), and the 2DEG is amenable to lo- 
cal probes including those which will allow the study of 
inhomogeneous states. In addition, the density and the 
polarity of carriers in graphene can be adjusted simply by 



changing the gate voltage; something that is not possible 
in conventional semiconductor heterojunctions. 

Here, we investigate the possibility of Wigner crystal- 
lization in graphene. Following the pioneering work of 
E. Wigner [1^, a Wigner crystal has become one of the 
extensively studied phases of conventional 2DEG at low 
densities. This state appears when electrons localize and 
form a crystal to minimize the potential energy, while 
paying the concommitant kinetic energy cost which arises 
from localization, as the density of carriers is lowered. 
Theoretical studies predict that Wigner crystallization in 
conventional 2DEG occurs at ~ 37, where r^QB is the 
mean interparticle spacing and as is the Bohr radius p^ . 
In the presence of a strong magnetic field, the kinetic en- 
ergy of electrons is quenched at filling factors z/ < 1; in 
other words, electrons can localize without paying the 
kinetic energy cost. Therefore, Wigner crystallization 
is facilitated by the magnetic field. The experimental 
study of the quantum hall effect in conventional 2DEG 
has shown the existence of the Wigner crystal phase at 
filling factors v < 1/5, deduced from vanishing Hall and 
longitudinal conductivity [l^ . 

The situation in case of graphene is not yet fully ex- 
plored. In the absence of a magnetic field, graphene 
shows metallic behavior over a large range of density for 
both polarities. In the quantum Hall regime it shows 
finite Hall and longitudinal conductivity, although only 
integer filling factors have been explored. Hence, it is 
appropriate to ask if Wigner crystallization occurs in 
graphene. In this note we show that Wigner crystal phase 
is not possible in graphene except in the presence of an 
external magnetic field. 

First, we present a heuristic argument for our result. 
Our discussion is limited to low temperatures T —^ 0; at 
finite temperatures, a translationally broken symmetry 
state will be destroyed by fluctuations in two dimensions 
[l5| . We know that the Wigner crystal phase is formed 
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as a result of the competition between the potential en- 
ergy Ep and kinetic energy Ek of the system. When 
Ep ^ Ei^ the crystallization occurs. We estimate the two 
energies in graphene as follows. Assuming localization of 
carriers at length scale I (which, in turn, is related to 
the density of carriers), the potential energy is given by 
Ep ~ {e'^/el)nA where n is the density of carriers, A is the 
area of the system, and e is the static dielectic constant 
of graphene. The kinetic energy, on the other hand, will 
be given by Ek ~ Tivci'^'^/ljnA where hvc = 5.8 eVA is 
the velocity of massless Dirac fermions. Hence, the ratio 
of these two energies is 



Ep 



h tvc 



(1) 



Two important observations thus follow immediately: i) 
The ratio is independent of the density of carriers, a situ- 
ation very different from that in the conventional 2DEG 
ii) The ratio depends on only two (tunable) material pa- 
rameters, e and vq. 

Thus, for graphene, we find that energetics of the sys- 
tem do not depend on the carrier density. In particu- 
lar, since it is known experimentally that graphene is 
not an insulator, it will not undergo Wigner crystal- 
lization to an insulator by reducing density of carriers. 
In the following, we present a microscopic calculation 
of the two energies, and obtain numerical estimates for 
suppressing the dielectric constant to increase the ratio 
Ep/Ek, which may drive the system towards a Wigner 
crystal. It is worthwhile to remember at this point that 
in conventional 2DEG, Wigner crystallization occurs for 
Ep/Ek = ~ 37. 

Now we will proceed with more detailed estimates of 
the kinetic and potential energy. In the ground state \F), 
the kinetic energy of the system is given by 



Ek = {F\Ho\F) = hvGj2k{F\al^ai,x\F) 



kA 



hvG 4:Ak% 



(2) 



where = 4 is the sum over the spin and valley degen- 
eracy, a]^_)^ creates an electron with momentum k (mea- 
sured from the Dirac point), and the Fermi momentum 
kp \^ related to the carrier density by n = k^jix. There- 
fore, the kinetic energy becomes 



^-^nvoAn^l\ 



(3) 



Similarly the potential energy is calculated as follows 

2 2_7r 
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(4) 



So their ratio is given by 

Ep 



h evG 



(5) 



and is independent of density of carriers as expected. For 
the experimental parameters, e = 4.7 and hvc = 5.8 
eVA, the ratio Ep/Ek — 0.32. Thus, kinetic energy is 
roughly three times bigger than the potential energy, and 
the carriers prefer to be delocalized in a uniform state. 
For Wigner crystallization, we need Ep/Ek ^ 1- 

Next we consider the dependence of this ratio on the 
two tunable parameters, e and vq- To this end, we first 
calculate the change in the static dielectric constant, 
e(q — > 0,w = 0) due to the spin-orbit coupling Aso in 
graphene [l^. The spin-orbit coupling in graphene is 
estimated to be as high as 4 meV 17]. Therefore, follow- 
ing the derivationn in ReflU, we calculate the dielectric 
constant for < Aso < 5 meV. 

At a; = 0, the dielectric functional e{q) can be written 

as 



e(q) = 1 + C/o(-7)n(q) 



(6) 



where UQ{q) — 2Tre'^/q is the two-dimensional Coulomb 
interaction in Fourier space and 11 (q) is the static bare 
particle-hole propagator 



n(q) 



-4 2^ l5k (q)l pA' _ pA , ■ (7) 

■ • ■ ^k+q -^k 
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Here, the factor of 4 takes into account of the spin and 
valley degeneracy, the index A = ±1 denotes the two 
bands, and the vertex factor is given by [lij . 



.9k'(q)l' = ;^(l + 



1 AA' sin Qfk+q sin at (fc + 9 cos ( 



|k + q| (8) 
-t-AA' cos Qfk+q cos Qfk), 



where 9 is the angle between k and q. is defined 
through, tanak = ^ 



Et-± 



A2 



h^v^k^ is the 



energy spectrum of the carriers in ± bands and f[E^] is 
the Fermi distribution function. 

At low temperatures T — s- 0, only the interband scat- 
tering contributes to e(q). Figure 1 shows the depen- 
dence of e on Aso evaluated numerically. We see that 
the dielectric constant decreases sharply with an increase 
in Aso- When the spin orbit interaction energy is zero, 
the dielectric constant is e = 4.7, whereas, e ~ 1.1 for 
Ago — 4 meV. The decrease in dielectric constant in- 
creases the energy ratio to Ep/Ek — 1-4. Even in the 
extreme case, e — 1 (no screening whatsoever) , the max- 
imum value of the ratio is Ep/Ek — 1-5, which is smaller 
than the ratio Ep/Ek = ~ 37 required for Wigner 
crystallization in a conventional 2DEG. 
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dergo Wigner crystallization as there is no concommitant 
increase in the kinetic energy cost associated with the lo- 
calization (This conclusion is robust for the n = Landau 
level, where the eigenf unctions are identical to those in 
the lowest Landau level of a conventional 2DEG). The 
critical value of the filling factor V(, at which the Wigner 
crystallization occurs will have to be determined by a mi- 
croscopic calculation, and is the subject of future work. 

This work was supported by DOE BES and LDRD at 
Los Alamos. We are grateful to T. Wehling, M. Katsnel- 
son, A. Lichtenstein, and R. Barnett for valuable discus- 
sion and advice. 



Figure 1: Static dielectric constant t{q ^ 0, = 0) as a 
function of the spin-orbit interaction energy Aso- We have 
used hvG = 5.8 eVA. We find that at Aso = 4 meV, e ~ 1.1 
approaches the minimum value possible e — 1. 

The second parameter vq can be tuned by putting 
strain on the system. In a tight-binding model, the veloc- 
ity VG is related to the tunneling amplitude t and the lat- 
tice constant a hy vg = 5ta/2. By straining the sample, 
the lattice constant increases linearly whereas the tun- 
neling amplitude decreases exponentially, thus effectively 
decreasing the velocity of massless carriers. How effec- 
tive this method will be in increasing the ratio Ep/Ek, 
depends upon the elastic properties of graphene; how- 
ever, it seems unlikely that the velocity vg can change 
by an order of magnitude that will be required for Wigner 
crystallization. 

Therefore, we conclude that Wigner crystal phase is 
absent in graphene. The kinetic energy is always com- 
parable to the potential energy, and both scale as n'^/^ 
with the carrier density. This peculiar dependence arises 
because of the linear dispersion of carriers. 

We note that these arguments will not be valid in 
the presence of a strong magnetic field. In that case, 
the system develops Landau levels with energies E,n = 
ihyG/lsV^m where Ib is the magnetic length. Each 
Landau level is macroscopically degenerate. Therefore at 
small filling factors <C 1, we expect the carriers to un- 
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